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Abstract 



The work presented in this thesis forms part of a program at the Paul Scherrer 
Institute (PSI) to construct a high intensity superthermal ultra-cold neutron (UCN) 
source based on solid deuterium (SD2) as UCN production medium. We carried out a 
set of experiments to gain a better understanding of the properties and the behaviour 
of solid deuterium as a cold neutron moderator and ultra-cold neutron converter. 

We present the measurements of the total neutron cross section as obtained by 
transmission studies with very cold neutrons and ultra-cold neutrons in solid deu- 
terium. The experimental set-up and the methods of data analysis are described and 
also the procedure of preparing the solid deuterium samples is given. The neutron 
transmission studies are supported by optical investigation of the crystal and by Ra- 
man spectroscopy. We have thus characterised the temperature dependence of the 
neutron transmission through solid deuterium and we have been able to identify the 
role that coherent neutron scattering plays for the investigated deuterium samples. 
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Chapter 1 
Introduction 



1.1 Physics of ultra-cold neutrons 

Ultra-cold neutrons (UCN) are neutrons of very low energy: the kinetic energy of a 
UCN is smaller than the effective Fermi potential Vp of a given material: 

where M„ = 939 MeV is the neutron mass, is the number density of atoms and 
h is the bound coherent scattering length of the material. Materials which are used 
in UCN physics are normally chosen to have a high Fermi potential, typically of the 
order of 10"'' eV . For example Be (Vp = 252 neV) and ^^Ni (Vp = 335 neV) are 
often incorporated as wall coatings in UCN equipment. UCN undergo total reflection 
at any angle of incidence from the surface of a material. As a consequence UCN can 
be trapped in bottles for times comparable with the free neutron lifetime of about 
887 s. Alternatively, UCN can be stored in magnetic traps: the UCN are confined 
to a trap region defined by a strong magnetic field that acts as a barrier through the 
interaction of the neutron magnetic moment and the magnetic field i2j . 

The fact that one can store and observe UCN for such long periods makes them 
an excellent tool to study fundamental properties of the neutron. Stored UCN offer 
the possibility of a precise measurement of the neutron lifetime that plays a key role 
for Big Bang nucleosynthesis and also is one of the main input parameters in tests 
of the Standard Model. Moreover, stored UCN are ideal for electric dipole moment 
(EDM) experiments. The existence of a permanent non-zero neutron EDM violates 
both parity (P) and time reversal (T) invariance and, through the conservation of 
the combined symmetry CPT, also violates CP symmetry (C refers to charge conju- 
gation). CP violation, which is believed to be responsible for the baryon asymmetry 
of the Universe, has been measured in the K° and the B'' mesonic systems. These re- 
sults are consistent with the description of CP violation in the Standard Model using 
the Cabibo-Kobayashi-Maskawa mass mixing matrix (CKM) ansatz. CP violation as 
described in the Standard Model is however at a level which is orders of magnitudes 
too low to generate the observed baryon asymmetry which suggests an existence of 
a mechanism for CP violation which is beyond the Standard Model. UCN can thus. 
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through a measurement of their EDM, form a very sensitive probe for physics be- 
yond the Standard ModeL Other fields where UCN can be used are also emerging: 
surface physics [3^ and the observation of quantum states in the Earth's gravitational 
field 14,^. 

Many of the experiments with UCN are limited by the statistics that is achieved in 
present sources. Next generation experiments that aim to improve these sensitivities 
rely on the development of new high intensity UCN sources. 



1.2 Methods of the UCN production 

Neutrons produced in a nuclear reaction in a spallation source or in a fission reaction 
in a nuclear reactor have energies of the order of MeV. These fast neutrons are mod- 
erated to thermal velocities by a moderator (e.g. heavy water) and can subsequently 
be slowed down into the cold regime upon interaction with a cold source (e.g. at the 
Institute Laue-Langevin (ILL) liquid deuterium at a temperature of 25 K acts as a 
cold source). To obtain neutrons with energy in the UCN region two methods can 
be used: 

• the conventional method which relies on extracting slow neutrons from the tail 
of the Maxwellian distribution of neutrons from the cold source, 

• a method in which cold neutrons are converted to UCN through a superthermal 
UCN production mechanism. 



1.2.1 Conventional UCN extraction 

The velocity spectrum of the neutrons in the cold source follows a Maxwellian velocity 
distribution. There will thus be a fraction, albeit a small fraction, of neutrons of UCN 
energy present in the low energy tail of the spectrum. The thermal spectrum of the 
moderator of temperature Tq is given by: 

p{v)dv = exp( -)— [1.2) 

a a'^ a'^ a 

where $o is the total thermal neutron flux and ot = >/2kBT})/M^ is the mean thermal 
velocity. Integration of ll.2l for the energies ranging from zero to the potential Vp that 
corresponds to the Fermi potential of the wall material gives the maximum UCN 
density: 

r ( W ( Vf \^ 

PUCN = / p{v)dv = -— -— 1.3 
Jo 3 a \kBToJ 

Ultra-cold neutrons can be extracted from the cold source through a vertical guide. 
Gravity is acting on the neutrons traveling upwards and decelerate faster neutrons 
to the lower energies. The mechanical deceleration can be used as well: a moving 
scatterer interacts with the neutrons and will carry away some of the neutron momen- 
tum. The number of UCN generated however is limited by Liouville's theorem: the 
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phase space density of the neutrons must be conserved if the particle move under the 
action of forces which are derivable from a potential. The final UCN density cannot 
be higher than the density at the cold neutron source because the compression in the 
neutron velocity space dilutes the neutron number density. 

The method described here is used at the UCN source operating at the Institute 
Laue-Langevin in Grenoble, France. It is in fact the strongest source of UCNs in the 
world. Extraction of the UCN from a liquid deuterium cold source is done through a 
vertical guide. Neutrons that travel the 17 meters upward through the guide and pass 
through a curved guide part at the end have kinetic energies in the very-cold neutron 
region. These VCN are guided into a turbine where they interact with moving mirror 
blades that slow them further down, by mechanical deceleration, to UCN energies. 
The details of the turbine can be found in Ref . jS] . The UCN density obtained in this 
set-up reaches 40 cm~^ in a storage bottle placed immediately outside the turbine. 
Further improvement in this density is not expected because of the limitation imposed 
by Liouville's theorem. 

1.2.2 Super thermal UCN production 

The principle of superthermal UCN production 

UCN can also be created from cold neutrons in a superthermal process. In this process 
neutrons interact with a medium that acts as a moderator to the cold neutrons and 
converts these neutrons into UCN by inelastic scattering. The medium must posses 
specific properties such as well defined energy levels and excitations that enable this 
down-scattering to take place. Superthermal UCN production was first proposed by 
Golub and Pendlenbury in a paper in which the UCN production in superfluid ^He 
through this mechanism was studied The understand the principle, consider a 
moderator with two energy levels, a ground state and an excited state, separated 
by an energy A. The moderator in the ground state can be excited by absorbing 
an energy A through inelastic interactions with neutrons of energy Eucn + The 
neutrons will thus be down-scattered to UCN energies. The inverse is also possible: 
the moderator falls back to its ground state by releasing an energy A and UCN 
can absorb this energy and thus be up-scattered out of the UCN energy region. The 
principle of detailed balance for this system leads to a relationship between the down- 
scattering rate (UCN production) and the up-scattering rate (UCN losses): 



where aup is the up-scattering rate of a neutron with energy Eucn and a down is the 
down-scattering rate of a neutron with energy Eucn + A. At low enough tempera- 
tures the up-scattering rate will be suppressed and the down-scattering will dominate 
resulting in a net production rate of UCN. It is noteworthy that this process is not 
limited by the Liouville's theorem for the neutrons alone. 

The UCN density that will be obtained with a system in equilibrium depends on 




(1.4) 
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the ratio of the production rate P and the loss rate 1/r: 

p = P-T (1.5) 

The loss factor is determined by a combination of loss mechanisms such as the neutron 
lifetime and losses specific to the design of the superthermal source such as nuclear 
absorption in the moderator, wall losses etc. 

Superthermal UCN sources 

There is a range of materials that could, in principle, be suitable as superthermal 
converters. Superthermal UCN sources based on superfluid helium and on solid 
D2 have already reached an advanced stage and have shown to be generating high 
densities of UCN. Furthermore converters such as CD4 |2| and O2 |H1 IH] are currently 
being studied and appear to be promising candidates for further progress in this field. 

The first superthermal converter that was proposed is superfluid helium [H]. One 
of the major advantages of superfluid helium is that the storage lifetime of UCN 
in such a medium is very long: ^He has got zero neutron absorption cross section 
and the up-scattering modes can be strongly suppressed by keeping helium at tem- 
perature of 0.5 K. UCN are produced in superfluid helium by down-scattering cold 
neutrons through creation of phonons. The dispersion curves of a free neutron and 
a phonon intersect, apart from at the origin, only at one energy of 1 meV. This 
means that, due to conservation of both energy and momentum, neutrons of only 
one wavelength, 9 A, can be down-scattered to UCN energies. The fact that the 
production of UCN can proceed through the interaction of cold neutron of one single 
energy only results in a helium UCN production rate about two orders of magni- 
tude lower as compared with solid deuterium where there is no such restriction on 
the incident neutron energy in the down-scattering process. Furthermore, efficient 
extraction of UCN from the superfluid has proved to be extremely difficult and thus 
most superthermal UCN sources are incorporated in single dedicated experiments (a 
neutron lifetime experiment at NIST jTU], neutron EDM experiments at LANL jlTj 
and ILL P2])- A superthermal UCN source based on superfluid helium that will feed 
several independent experiments is being studied only at Osaka P^ . 

After the first work on ^He UCN sources, other materials emerged as a possible 
medium to create UCN. Most noticeable is the development of the solid D2 based 
sources (that have an absorption lifetime of 150 ms, which is about four orders of 
magnitude lower than that for He) together with the idea of the "thin film superther- 
mal sources" ^^Ej. The thin film superthermal source is based on the notion that 
the UCN density in a storage volume that contains converter material is independent 
of the actual volume of the converter material. The use of the pulsed neutron source 
in combination with the converter with a short neutron lifetime was first pointed out 
by Pokotilovski P^l- In this scheme the moderator is connected to a UCN storage 
volume only for a short time. After the neutron pulse has occurred a valve closes 
separating the conversion and storage regions, this enables the UCN to survive for 
long times as they are no longer interacting with the absorbing converter material. 
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Superthermal UCN sources based on solid deuterium are currently under intense 
investigation at PSI [HI and at Los Alamos [HI US 1201 IHI • The production of UCN 
using sD2 has first been studied at PNPI Gatchina where a gain factor of 10 in UCN 
production was obtained compared to the UCN production using liquid deuterium 
source [221 12n] • At Los Alamos a prototype SD2 based UCN source was operated and 
UCN densities in excess of 100 cm~'^ have been reported. The UCN source being built 
at PSI, also based on SD2, aims at achieving densities of 3000 cm~^ in a 2 m^ volume. 
A detailed discussion about SD2 as a UCN moderator is presented in the next chapter. 



Chapter 2 

sD2 as a super-thermal source 



A good superthermal converter needs to possess a number of properties that are 
directly related to the creation and losses of UCN. In this chapter the properties of 
deuterium are described with the focus on its use as superthermal UCN converter 
material. Some of the properties of deuterium are already quite favourable: deuterium 
has a relatively low neutron absorption cross section (the absorption cross section for 
thermal neutrons is 519 |J.b) and it is a light element. Also, deuterium has a large 
down-scattering cross section because the phonon density of states has a large overlap 
with the cold neutron spectrum. 



2.1 Neutron scattering in SD2 
2.1.1 Coherent and incoherent scattering 

To investigate the suitability of using solid deuterium as a superthermal UCN con- 
verter it is important to understand the neutron scattering from this material. The 
interaction potential between a neutron and a rigid array of N nuclei in the first Born 
approximation is j24j : 

nr) = ^fj2^i6{r-Ri) (2.1) 

where Mn is the neutron mass, bi is the bound nucleus scattering length for the /th 
nucleus and is the position vector of the ith nucleus. The final state ipii' and the 
initial state ip^ of a neutron, defined by the wave vectors k' and k, respectively, are 
given by plane waves. For a large box of volume the states ip^ and ip^r are: 

exp(ik-r) . , ^ h'^k'^ 

^k= with energy E = -— (2.2) 

and 

exp(zk'-r) , _ h'^k''^ . . 

^k' = -J with energy E = — — (2.3) 

La /iWji 

The initial state of the system composed of the incident neutron and the scatterer is 
described by: 

\k) |A) = \kX) (2.4) 
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where |A) denotes the state of the scatterer with energy Ex. 

The transition rate from the initial state ■i/'k to the final state V^k' is given by the 
Fermi's golden rule 



27r 

dWk^k' = -- 
n 



p^,{E') 6{E + Ex-E' ~Ey) (2.5) 



where V refers to the interaction potential that causes the transition. py.,[E')dE' is 
the density of the neutron final states per unit energy interval 

pAE') = ^^^dn. (2.6) 

The cross section is obtained by dividing the transition rate dW by the incident flux of 
neutrons {hk/{MnL^)). Using the Dirac notation we can write the double differential 
cross section as: 

"^'"^ ^'\{k'X'\V \kX)\^ 5{hiu + Ex-Ey) (2.7) 



dQdE' k 
where hu is the defined by: 

hu = E - E' = -—(e - k'^) (2.8) 

Inserting expression 12.11 for the scattering potential one obtain the differential cross 
section per solid angle per atom: 

= M E^^K^'I E ^^P(^Q ■ ^iM'^i^^ + Ey) (2.9) 

A A / 

where = (k — k')h is the momentum transfer at the scattering, px repre- 
sents the weight for the state |A), which is proportional to the Boltzmann factor 
exp{-Ex/kBT). 
We can write equation 12.91 as: 

7 = E WA| exp(-^Q ■ R,)|A') 



dQdE' N k , 

A A l.l 



x(A'|exp(iQ-R,)|A) 6{huj + Ex - Ey) (2.10) 



The quantity h^ihi is the average of 6^*6; over a random nuclear spin orientation and a 
random isotope distribution and is described by: 



hp, = \w+5,Ar\-\W) (2.11) 

Substituting 12 . 1 ll into 12 . 1 01 the cross section can be split in two parts: 
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where the coherent cross section is: 

\ / coh X'X I 

^jblhScoh{Q,u;) (2.13) 

and the incoherent cross section 

V / mc X'X I 

^|6L'5mc(Q,^) (2.14) 

S{Q,u!) represents the response function or dynamic structure factor and huj is the 
energy lost by the neutron in the scattering process. The response function is given 
by the Fourier transform of the density correlation function G{r, t) which represents 
the probability that, given a particle located at the origin at time t — 0, any particle 
is found at position r at time t. 

Elastic scattering 

If the scattering system has no internal structure (no excited states) then the energy 
of the scattered neutron is identical to that of the incident neutron. Such scattering 
is referred to as elastic scattering. For elastic scattering Ex — Ex' and the total 
coherent and incoherent cross section reduces to: 

a^h = 47r6^„^ (2.15) 

ainc = 47r6L (2-16) 
In the coherent scattering there is a strong interference between the waves scattered 
from neighbouring nuclei. Thus the coherent scattering is present only if strict geo- 
metrical condition are satisfied. In the incoherent scattering there is no interference 
at all and the cross section is isotropic. 

When a neutron with spin 1 — 1/2 interacts with a nucleus of a deuterium atom 
of spin 7 = 1 it can do so in states of total spin S — 3/2 or S — 1/2. Therefore the 
scattering lengths fe*^"*"^ and 6'-"^ should be associated with these two states. We have 
4 substates of spin 3/2 and 2 substates of spin 1/2, so the statistical weight of the 
8=3/2 state is 2/3 and that of 5 = 1/2 state is 1/3. The scattering lengths bcoh 
and bine are thus averaged over these two scattering modes: 

bcoh = lb(^^ + lb(-^ (2.17) 

and 

knc = ]Jl\b(+^-bcoh\^ + ^\b(-)-bcoh\' (2.18) 

Using the measured values 6*^+^ = 9.5 fm and b^~^ = 1.0 fm we find acoh = 47r6^o^ = 
5.59 b and ainc = 47r6i = 2.04 b. 
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Scattering from static inhomogeneities 

If the neutron scattering takes place on nuclei that move we are dealing with an energy 
dependent function S{Q,lj). For elastic scattering with a static target, where hu is 
zero, the response function S{Q,uj) = S{Q), where S{Q) is called "structure factor". 
All the information about the static distribution of atoms in the scattering system is 
then represented by the structure factor. The scattering from static inhomogeneities 
with radius R and uniform density difference to the surrounding medium Sp = p — Po 
leads to |n|: 

S{Q) = {4nR'6pr (^^^y (2.19) 

where ji{x) is the spherical Bessel function of order 1. Thus the total coherent elastic 
cross section is 



■'-"i-H.^U^S.f r Wml,Q (2.20) 



where 61 is the detector aperture. The integral in 12.201 is 



r bi{QR)]\ Q ^ g^^j^Q^^ ^ hjiikRe,) + jlikRe,) - jl{2kR) - jl{2kR)] (2.21) 

In section 14.2.11 this model will be used to interpret the transmission data from the 
VCN neutron experiment. 

2.1.2 Inelastic neutron scattering from deuterium molecules 
Deuterium molecule 

The deuterium molecule consists of two deuterons of spin 1 and thus can form states 
of total nuclear spin S = 0, 1 and 2. State of S = and 2 refer to ortho-deuterium 
and state of S = 1 refers to para-deuterium.^ The molecular dynamics of deuterium 
is determined by the motions which are represented as the translation of the center 
of mass, the inter-molecular oscillations and rotations. To excite the first vibra- 
tional state (vibrational quantum number n = 1) of deuterium molecule an energy of 
386 meV is required and thus for the low energy neutrons^ the vibrational excitations 
do not contribute to the scattering cross section. The rotational energy spectrum for 
diatomic molecules with the moment of interia / = mR^ is given by: 

Ej = ^^J{J + l) (2.22) 



""^At room temperature the deuterium is a mixture of 66.7% ortho-D2 and 33.3% para-D2. Higher 
ortho-D2 concentration up to 98.5% can be achieved by using the ortho-para converters |27j . 
■^Cold neutrons have energies below 10 meV 
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Using for the deuterium molecule the values m = 3.75 GeV and R = 0.37 A ^ the 
energy becomes Ej = 3.75 J( J + 1) meV. 

The state with total nuclear spin S = 0,2 can only have even values of rotational 
quantum number J, and for the para state {S = 1) only odd values of J are allowed. 
The population of a state, Nj, with the rotation quantum number J follows the 
Boltzmann distribution: 

^ = (2J-H)exp(-A,) (2.23) 

where T is the temperature of the system. At temperatures where deuterium is in 
the solid phase T < 18 K the molecules are in their lowest rotational state, which is 
J = for ortho-deuterium and J = 1 for para-deuterium. The transitions between 
even and odd J states are prohibited by the spin selection rules in the case of a non- 
interacting gas. For the solid phase such a transition is possible due to a spontaneous 
spin flipping, however, such a conversion is extremely slow compared with the time 
scale of experiment [HI EH] ■ 



Properties of the SD2 lattice 

To evaluate the cross section for solid deuterium one has to consider the dynamics of 
the lattice. Many properties of solids can be successfully described by a simple con- 
tinuum Debye model. The dynamical properties of the solid lattice are characterised 
by the Debye temperature 6d, which defines an effective maximum phonon frequency 
uj£) in the solid: 

en = ^ (2.24) 



For solid deuterium the Debye temperature is 6'/) = 110 K 
The normalized phonon density of states according to the Debye model is: 

Z(uj) = ^ (2.25) 

The position vector of the /th atom in the crystal of N atoms can be written as: 

Ri = pi+ui{t) (2.26) 

where pi represents the equilibrium position of the Zth atom. The displacement ui{t) 
of the atoms from their equilibrium configuration for the harmonic lattice can be 
written as ^: 

Mt) = E [^«(^' "iK^ + q)<J (2.27) 

s,q 

where q and the phonon annihilation and creation operators 



,) = l^s,q-l) (2-28) 



^The separation between the two nuclei in D2 is 0.74 A 
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and 



a7Jn,,„)= a/^^^^TT l^s.q + 1) (2.29) 



= ^ImMJM)^'^"'^ exp(^q.p) (2.30) 

u!s{(i) represents the normal- mode frequencies, the index s=l,2,3 denotes the three 
solution for uj for each q [21] and 7s (q) is related to the polarization vector. We 
assume that all the atoms have the same mass M. Substituting 12.261 and 12.271 into 
the matrix element of 12.141 gives 

(A'|exp(iQ-Ri)|A) =exp(iQ-pi) 

X 1[{X'\ exp {^Q • [e,,(^, q)«.,q + q)<J } |A) (2.31) 

s,q 

Applying a Taylor expansion of the exponential function in the operator we get 
exp {iQ • [l^i{s, q)a,^q + $*{s, q)a+J } = 1 + iQ ■ [Ci{s, q)a,^q + ^*(s, q)a+q] + 

'l [(Q -eOXq + (Q -^'O'Kq)' + IQ ■6.|'(as,q<q + <q«^,q] + " " " 

(2.32) 

Since q and a+q annihilate and create, respectively, phonons, we see that the second 
term changes the phonon number by one while the third term changes the phonon 
number by two or zero. 

To calculate the one-phonon contribution to the neutron cross section we can use 
UTM ?rM and ran and write 1^11 as: 



(A'l exp(zQ ■ R,) I A) = exp [-W{Q)] ^ / ^^^^^^^^ " 7(^, q) 

exp [ipi ■ (Q — q)] y^ris^q + 1 phonon creation 
X { ' (2.33) 

exp [ipi ■ (Q + q)]^7^s,q phonon annihilation 

where the term exp[— W^(Q)], the Debye- Waller factor, describes the zero phonon 
expansion. For a cubic lattice Vr(Q) is 24^: 

WiQ) = ^Q2(u2) (2.34) 
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where (u^) is the mean squared displacement of a nucleus. 

Using the incoherent double differential cross section of 12 .141 becomes: 

{[^^(q) + 1] 6{hLJs,q — huj) phonon creation 
(2.35) 
ns(q) 5{hiJs,ci + hcu) phonon annihilation 

Finally, after replacing the summation over s,q by the normalized density of states (25 
one obtains: 

^ ^ ^ ''LTexp|-2W'(Q)]{|q-7,(q)P}.„ 



dVtdE' L.^ 2M '""^k 



mc 



ryf \ ( niuo) + 1 if C(j > 0, phonon creation 
Z(uj] \ ^ ' — 1 f 



x^^<^ (2.36) 
^ n{uj) if cij < 0, phonon annihilation 

where hu is the energy lost by the neutron. The subscript av in 12.361 stands for the 
average over a surface in Q with constant u and in cubic symmetry this average 
is equal to ^Q^- The coherent scattering is usually small because it is restricted 
by the energy and momentum conservation and very often it is sufficient to use 
the "incoherent approximation" , which consists of using 12.361 with 6?^^ replaced by 
Koh + ^Inc- The vahdity of "incoherent approximation" is discussed in 0. 

We can now apply the resulting formula l2.36l to the solid ortho-deuterium system 
and find an expression for the cross section. Deuterium crystallises in the hep struc- 
ture with an intermolecular separation of 3.7 A. At low temperatures ( T < 18 K ) the 
vibrational n and rotational J quantum numbers are zero. Thus the ortho-deuterium 
molecule is spherically symmetric and may be considered as a single particle with a 
neutron scattering length 2bcohjo{^) 30j where bcoh is the coherent scattering length 
of the deuteron (Eg. 12.171 ). jo is a spherical Bessel function, a denotes the separation 
between the deuterons in the molecule (a = 0.74 A) and Q is the momentum transfer: 

Q = y/k^ + {k'y -2kk' cose (2.37) 



where 6 is the angle between k and k'. Following Liu et al. [2H] we can use a cubic 
lattice to describe the solid deuterium structure so that 



/ 79 \ 1 phonon r / ^ \ n 2 7 / 



D2 



^bcohJo [ -j- 



^exp[-2H^(Q)] 



„, X f niio) + 1 if cij > 0, phonon creation 
x^<! (2.38) 
n{uS) if C(j < 0, phonon annihilation 



The population number niiS) of phonons with the frequency 



a; = -^Q^ (2.39) 
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is given by Bose-Einstein statistics 



1 



niu] 



- 1 



(2.40) 



Figure 12.11 shows the up-scattering cross section for VCN calculated from 12.381 
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Figure 2.1: The calculated VCN up-scattering cross section vs. temperature of solid 
D2 with an initial neutron energy of 7.9 ■ 10~^ meV, corresponding to a velocity of 
66 m/s or a wavelength of 60 A. 



Chapter 3 
Experiment 



3.1 Experimental setup 

The experiment was carried out at the Intitute Laue-Langevin (Grenoble, France) on 
the VCN and UCN beam-hne during 2003. A total of 2 reactor cycles (reactor cycle 
135 and 136, 50 days each) had been allocated for these experiments. Reactor cycle 
135 was devoted to the VCN measurements and reactor cycle 136 was devoted to the 
measurement on the UCN test position. 

3.1.1 General overview of the VCN and UCN setup 

Figure 13.11 shows a schematic view of the experimental setup at the VCN beam 
position. Neutrons with velocities of 20 m/s to 200 m/s pass the setup as discrete 
packets created by a neutron chopper and are detected in a position sensitive micro- 
strip gas counter. The SD2 crystal was grown in the target cell with a volume of 
40 cm^. The cryostat with the target cell was placed at the end of the neutron flight 
path between the chopper and detector. During the experiment the detector was 
located at two positions: "front" and "extended" . For the "front" detector position 
the total length of the flight path^ is 1434 mm. The collimator assembly of a length of 
1260 mm consists of: the 20 mm diameter boron collimator placed directly after the 
chopper, the 20 mm diameter cadmium collimator placed before the target cell and 
four Cd/Cu collimators (Fig. 13. 2j) of total length 150 mm placed 600 mm after the 
chopper. The reasons for having the detector at two position are the calibration of 
time-of-fiight spectrum and the scattering measurements at small scattering angles. 

A similar setup was used for the UCN experiment (Fig. 13. 3|) . For these mea- 
surements the cryostat with the SD2 crystal was placed in the flight path of length 
1247 mm directly after the chopper. For the neutron scattering measurements a 
collimator with a length of 30 mm and a diameter of 20 mm was placed upstream 
of the chopper. After passing through the cell neutrons enter a neutron guide that 
consist of a movable tube. With this tube setup, scattering measurements at various 
acceptance sohd angles could be performed. 



Distance between the chopper and the middle of the neutron detector 




Figure 3.1: The scheme of the VCN setup, all dimensions are given in mm. The parts 
of collimation assembly are marked with red. The distance between the middle of 
the cell and the detector is 110 mm. 
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CiJ06O 




Neutron beam 



CU0 2O 



Cd 30 

Figure 3.2: The collimator stack which consists of four Cd/Cu colhmators, all dimen- 
sions are given in mm. The given diameters are the inner diameters. 
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3.1.2 The cryostat and the target cell 

In this section the general concept of the target cell and the cryostat will be described 
together with the method of preparing the crystal. Further details can be found 
elsewhere ^7\. Figure EiH shows the cryostat designed to obtain temperatures down 
to 5 K in the SD2 cell. These temperatures are required to study solid deuterium 
transmission properties for the superthermal UCN source. 




Figure 3.4: The SD2 cell mounted on the cryostat 



The cell is made of an aluminum alloy AlMgs and has four windows: two optical 
windows and two windows for neutrons. The optical windows, made from sapphire, 
serve for Raman spectroscopy as well as for optical investigation of the crystal. The 
neutron windows, made from aluminium, are 0.15 mm thick and have 26 mm diam- 
eter. The distance between the entrance window and the exit window is 10 mm. 
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Method of freezing 

The cell is designed in such a way that the crystal is frozen from the liquid phase. A 
high-resolution digital camera permanently monitors this process. Gaseous D2 kept 
in the aluminium storage container is passed through the gas system into the cryostat 
and is subsequently liquefied in the target cell. The process of freezing is carried out 
at a temperature close to the triple point (18.7 K), this procedure takes about 12 h. 
The crystal that is obtained under these conditions is fully transparent and clear. In 
Fig. 13.51 a series of pictures shows how the crystal is being grown from the melt. 




Figure 3.5: A series of pictures showing the crystal being grown from the melt 



The speed at which the deuterium is frozen is a very important parameter in this 
process. Figure ITT)! shows crystals grown at two cooling speeds and demonstrates the 
dependency of the crystal quality on the speed of cooling. Because of a temperature 
gradient that arises during the filling of the cell the deuterium in the upper part of the 
cell will be frozen before the complete crystal is formed. This will affect the flow of the 
rest of the deuterium as it moves into the cell and can cause problems to the formation 
of the remaining part of the sample. However, we have developed an annealing 
technique in which the deuterium is kept at a constant temperature for several hours. 
This annealing procedure has greatly improved the crystal homogeneity. Figure 13.71 
shows the annealing process of a crystal as it is grown in the cell. Once the crystal 
reached a temperature of 18 K it was cooled down further to a temperature of 5 K. 
This last step in the process of creating the flnal form of the solid deuterium was 
performed at a reduced cooling speed: in most cases with a typical cooling speed of 
2 K/h. 

The interpretation of the data is based on the measured neutron transmission 
through the sample. It is thus important to know accurately the amount of deuterium 
of the neutron beam passage. The deuterium concentration as a function of position 
in the cell is obtained by analysing the pictures as recorded by the digital camera 
(section Emi). 
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Figure 3.6: Two crystals at 5 K that have been grown at two coohng speeds. The 
crystal pictured on the right was grown at a speed of 2 K/h and the crystal pictured 
on the left was grown at a speed of 14 K/h. 




Figure 3.7: The annealing process of a crystal, the white line marks how the 'empty 
part' of the growing crystal is filled with time. Shown is, from left to right, 1 - start 
of annealing, 2 - after 5 h, 3 - after 15 h, 4- after 30 h, 5 - stop of annealing (40 h) 



3.1.3 The Raman spectroscopy setup 

The system which is used for Raman spectroscopy consists of three main parts: the 
argon-ion-laser, the Raman head, and the Raman spectrometer. A detailed descrip- 
tion of the setup and method of the data analysis can be found in Ref. |27j . 

Raman spectroscopy has been used to investigate the ratio of ortho-D2 to para-D2 
as well as the crystal structure. Hydrogen contaminations in form of HD or H2 can 
also be measured. The ortho-para concentration can be calculated from the intensity 
ratio of the So(0) and So(l) lines of the rotational Raman spectrum f section l^!T|) while 
the lattice structure of the solid deuterium can be determined from the So(0) signal 
which produces a multiplet unique to the crystal structure. From the relative line 
intensity within the multiplet one can deduce the crystal (crystallites) orientation 
with respect to the direction of the Raman collector f section 14.11) . Figure ITHl shows 
the Raman head mounted on the cell. 
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Figure 3.8: The Raman head consists of the laser head (upper) for the hght input 
and Raman collector (lower) for the scattered light. The Ar laser light is coupled to 
the laser head via optical fiber. 



3.1.4 The Camera 

Throughout the measurements the digital camera is permanently mounted on the cell 
(Fig. 13. 9|) . The camera allowed to monitor both the optical properties and the total 
amount of deuterium that is formed. More information can be found in Ref. ;27j . 

The pictures generated by the camera are used to study the optical difference for 
a specific crystal at different temperatures and prepared in different ways. 

3.1.5 VCN/UCN detector and DAQ 

The same neutron detector and DAQ system has been used for the VCN transmission 
and the UCN transmission experiments. Although the efficiency of the detector varies 
with the neutron energy, relative measurements are obtained directly by comparing 
the neutron transmission through the different samples and measurements taken with 
the cell being empty. 

Detector 

The neutron detector, a so-called "BidimSO" developed at the ILL, is a 2-dimensional 
micro-strip detector that consist of an MSGC plate placed in a sealed cylindrical metal 
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container filled with a mixture of 50 mbar of ^He and 1 bar of CF4. The entrance 
window to the detector is a 100 |J.m thick Al foil. Because of the 50 mbar overpressure 
the detector window was slightly bulged but this has no noticeable effect on the 
final results. A part of the detector is shown in Fig. 13.101 the MSGS plate and 
an aluminium diaphragm are mounted on a stainless steel flange that is equipped 
with a valve through which the gas can be filled. The two electrical feed-throughs 
mounted on a flange are used to apply the voltage onto the electrodes and to extract 
the neutron signal. 

The MSGC is made of a semiconductive glass plate with a chromium pattern of 
strips, the total sensitive area is 80x80 mm^. The 10 |^m wide anodes are engraved on 
the front side of the plate, the 980 (xm cathodes are placed on the rear side of the plate 
orientated perpendicular to the anodes. The spatial resolution for neutron detection 
is about 1.5 mm. The X and Y position of the detected neutron is determined from 
pulse height asymmetries of signals from the electrodes in the following way: 

X = ^^.Acm (3.1) 

y = ~ ■ 4 cm (3.2) 

A2 + AS 

The position of the signal on the detector was obtained by performing measure- 
ments with a cadmium mask (Fig. 13. lip placed on the front flange of the detector. 
Figure 13.121 shows a two dimensional picture from the detector with the Cd mask 
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mounted. The distance between the window and the MSGC plate is 30 mm (for a 
flat window). 

The high voltage applied between the entrance window and the anodes and cath- 
odes is 1200 V. Neutron are detected through the reaction '^He(n,p)^H. The energy 
release in this reaction is 764 keV and the cross section is 5330 b for thermal neu- 
trons. Scaling the cross section to VCN and UCN energies gives a reaction cross 
section of the order of mega-barns and forms a very efficient basis for neutron de- 
tection. The proton and triton emerging from the reaction will create electrons by 
ionising the gas mixture. The electrons drift in the electric field toward the MSGC 
plate. On approaching the glass plate these electrons are accelerated in very strong 
electric field in the vicinity of the anodes leading to secondary ionisation and creation 
of avalanche. The electrons penetrate the glass plate and produce electron-hole pairs 
along its track, the number being proportional to the energy loss. The electric field 
applied between the anodes and cathodes separates the pairs before they recombine; 
electrons drift towards the anode, holes to the cathodes on the rear side. Due to 
the semiconductive properties of the glass, the signal from the rear side is equal to 
the one on the front side. The electric pulses thus generated are passed through a 
pre-amplifier before being converted by ADCs^. A typical pulse height spectrum is 
shown in Fig. 13.131 




Figure 3.10: The MSGC plate and an aluminium diaphragm of the "BidimSO" de- 
tector. 



2port Modules and the MPA-3 BASE module 
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Figure 3.11: The cadmium mask that is used to determine the position of the signal 
on the detector. 
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Figure 3.12: A two dimensional picture from the detector with the Cd mask mounted, 
the numbering is the same as for Fig. 13.111 
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Figure 3.13: Pulse height spectra as obtained from the measurement with the empty 
cell at the "front" detector position. 
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DAQ and data format 

An MPA-3 Multiparameter System was used for data acquisition. The complete 
description of this system can be found in Ref. The main hardware modules of 
MPA-3 are: 

• An MPA-3 PCI Card with an internal 20 MHz clock reset by a TTL signal from 
the chopper. 

• An MPA-3 BASE module with four ADC ports for the connection of external 
ADCs and two BNC type auxiliary connectors AUX 1 and 2. One of this is 
used to accept a TTL signals from the chopper. 

• ADC Port Modules: an interface for the 4 ADCs. 

Data is written in a list mode file format and stored on the RAID disk system. A 
detailed description can be found in Ref. jSI]. Each event contains the information 
on the ADCs channel, 16 bit ADC data and real time clock RTC value in three 16 
bit words rtcO, rtcl and rtc2. The RTC value is calculated according to the formula 
(rtc2*65536-|-rtcl)*65536+rtc0. It starts from a preset value which was rtc0=65528, 
rtcl=65536, rtc2=65536 and counts down in 50 ns intervals. 

3.1.6 The neutron chopper 

All the measurements of the neutron transmission through the crystal have been 
performed at specific neutron wavelengths selected using a neutron chopper and time- 
of-fiight techniques. For the UCN and VCN measurements different choppers have 
been used. 

The VCN chopper 

The VCN chopper uses a rotating disk made of aluminium coated with a neutron 
absorbing Gd and ^LiF layer and it has two slits through which neutrons can pass. A 
20 mm diameter boron collimator is placed directly after the chopper and is centered 
on the beam axis. The chopper was operated at a rotation frequency of z/ = 7 Hz. 
A 5 fxs long TTL signal is generated before the slit reaches the neutron beam center. 
This signal triggers the data acquisition system. The cross section area of the beam 
that is cut by the chopper slit and the collimator at a given time defines the chopper 
opening function. This can be calculated by integrating the overlapping area of the 
slit and collimator as the slit moves through the beam area. The width of the TOF 
channel (80 yis) defines the time step for this calculation. The chopper opening 
functions have been calculated for the nominal chopper frequency of = 7 Hz and 
also for a chopper frequency of z/ = 3 Hz at which some measurement were carried out. 
These chopper functions are plotted in Fig. 13.141 The transmission profiles obtained 
by the opening functions can play an important role in the interpretation of the TOF 
spectra that are measured: the transmission profiles have been used to deconvolve 
the raw data. The details of the deconvolution are given in Appendix IA.1.21 
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Figure 3.14: The VCN chopper opening profiles obtained for two chopper frequencies: 
3 Hz and 7 Hz. 



The UCN chopper 

The UCN chopper consists of three rotating disks made from polycarbonate: one 
disk rotates slowly while two other disks rotate at the higher speed. The ratio of the 
speeds of the disks is fixed at 1:6. The 'slow' disk opens a time window of 22 ms with 
a typical operation frequency of about 3 Hz to 4 Hz and produces a neutron pulse. 
The two fast disks sharpen the rising and falling slope of the pulse. The disks in the 
chopper are driven by a stepping motor, the stepping motor is computer-controlled 
through a "PCI 7324" card and Lab View software. The trigger signal needed for the 
TOF measurements is provided by the light being reflected by one of the reflection 
sensors and being detected by infrared diodes. The reflector is a 7 mm long piece of 
aluminium mounted on the wheel of the slow disk, which together with the rotation 
speed determines the length of the pulse. To correct for this effect, a shaping stage 
was applied. The length of the TTL signal, after passing this stage is 5 |xs. A full 
chopper turn takes 864 steps. The chopper was routinely operated at z/ = 3 Hz. The 
stepping motor can be run at a slow enough speed that allows for an experimental 
determination of the chopper opening function. Fig. 13.151 shows the neutron yield 
(and thus the chopper opening) as a function of time (blue points). The results with 
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the best statistic are obtained at a chopper speed of 2 steps/s which corresponds 
to 2 mHz. The curve thus measured has been fitted with a polynomial function, 
Fig. 13.1 51 shows the resulting curves. In order to get the chopper opening function for 
the typically used speed of 2500 steps/s, we used the obtained polynomial parameters 
to evaluate appropriate values and the result of this procedure is presented in the 
Fig. 13.161 The effect of convolution of the obtained opening function with the raw 
UCN TOF spectrum was calculated to be negligible in the analysis of the data. The 
stepping motor can also be operated in a way so as to obtain a precise calibration of 
the TOF spectra. This technique is described in Appendix IA.2.21 
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Figure 3.15: The chopper opening function for the chopper speed 2 steps/s (2 mHz) 
plotted together with a fitted polynomial function (red line). The bin width is 100 ms. 
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Figure 3.16: The chopper opening profile for a chopper speed of 2500 steps/s (3 Hz). 



Chapter 4 

Results and data analysis 



A typical set of data consist of neutron transmission measurements through hquid 
and sohd deuterium. It took from 4 to 7 days. Measurements with the deuterium 
in the gaseous phase were also carried out in some cases. Once the transmission 
measurements of the deuterium in the liquid phase had been completed, the crystal 
was grown from the melt and was subsequently annealed f section l3.1.2j) . The trans- 
mission through solid deuterium is measured as a function of the crystal temperature 
range from 5 Kto 18 K in 2 K steps in most cases, from 18 K. Simultaneously Ra- 
man Spectroscopy was performed as well as an optical investigation with the digital 
camera. For the "front" detector position, two sets of measurements with deuterium 
containing a high ortho-D2 concentration were carried out and one set of measurement 
with normal deuterium was done. We also performed a few additional measurements 
testing in what way the method of freezing influences the neutron transmission. This 
section focuses on regular measurements. Others are only shown to emphasise some 
important effects that have been observed. 

The analysis of one deuterium crystal consists of the following: 

• Determination of ortho-D2 concentration from the Raman spectra (usually from 
the liquid measurements). 

• Optical inspection of the measured crystal using pictures recorded by the digital 
camera. 

• Analysis of the So(0) line splitting in the ortho solid which gives information 
about the relative orientation of the crystal (or crystallites). 

• Analysis of the neutron data including measurements with liquid, solid at var- 
ious temperatures and an empty cell. 

4.1 Deuterium properties 

To investigate the ratio cq of ortho-D2 to para-D2 deuterium in the measured sam- 
ple the Raman spectroscopy method has been used. Figure 14.11 shows spectrum of 
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rotational excitations for liquid ortho-deuterium with cq = (98.7 ± 0.2)%. The ortho- 
deuterium concentration is determined by comparing the Raman line intensities for 
the J=0— >-2 and J=l— >-3 transition (lines So(0) and So(l) respectively) and taking 
into account the proper transition matrix element. The Ar II (493.3209 nm) line 
which is seen in Fig. 14. H is produced by the laser. 





80000 






1 1 1 1 1 1 1 


1 ' 1 


' 1 




70000 




s„(0) 

1 








60000 












E 
o 


50000 












/(0.5 


40000 












Counts 


30000 
20000 




j 


1 Aril 

\ i 




s„(i)^ 




10000 






— — -_-_--/V______. 




i ■ 






1 




1 . 1 ■■ 





160 180 200 220 240 260 280 300 
Raman shift [ cm ^ ] 

Figure 4.1: A Raman spectrum of liquid 0-D2 ((98.7 ± 0.2)%). The spectrum contains 
data from 14 single files with 200 s exposure time each. 

Figure 14.21 shows the Raman signal for solid ortho-deuterium at temperature of 
18 K. The So(0) line is split into a triplet which is a signature of the hep structure of 
the crystal The hep symmetry of deuterium crystal makes the J = 2 level form 
three energy bands belonging to m = ± 1 (a) m = ± 2 (/?) and m = (7). Change 
of the So line shape with crystal temperature is shown in Fig. 14.31 

The pictures which were taken with the digital camera allow us to study the optical 
transparency of the deuterium crystal at various temperatures. It appears that only 
carefully prepared crystal at temperature of 18 K is fully transparent (Fig. 14. 4j) . 
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Figure 4.2: A Raman spectrum of solid ortho-deuterium at 18K. The spectrum con- 
tains data from 14 single files with 200s exposure time each. 




Figure 4.3: The J=0^2 transition for the deuterium crystal at various temperatures. 
The relative m^ = 0, 1, 2 intensity ratio is changing with the temperature which 
indicates that the crystal orientation is changing as well |32]- 



34 



CHAPTER 4. RESULTS AND DATA ANALYSIS 



4.2 Neutron data analysis 

The neutron total cross section was evaluated from the neutron beam attenuation 
in the deuterium sample. The neutron count-rate A^o has been measured for the 
beam passing through an empty target cell and the neutron count-rate A^i for the 
beam passing through the target cell filled with deuterium has been measured. These 
count-rates are related through: 

Ni = Noexpi-patotl) (4.1) 

where p is the density of the sample, / is the length of the cell and atot is the total cross 
section. The density of a deuterium sample was calculated from the molar volume 
Vl, taken from Ref. [221; in the following way: 

(4.2) 

where is Avogadro's number. The molar volume Vl is equal to 23.21 cm^/mol for 
liquid deuterium at 19 K and 20.34 cm^/mol for solid deuterium. 

A single run takes normally 1800 s. In most cases data were taken for 2 hours 
at one temperature. The empty cell transmission measurements have been also per- 
formed in the range from 5 K to 20 K. No noticeable temperature dependency has 
been found. 



4.2.1 VCN data 

Figure 14.51 shows neutron total cross section versus neutron wave vector k that have 
been measured with liquid deuterium and with solid deuterium at 18 K and 5 K. The 
solid deuterium data were obtained with the same crystal with a high ortho-deuterium 
concentration of (98.7 ± 0.2)%. 

The fitting was done with the MINUIT package following two minimisation 
methods: Simplex ^3] and Migrad [SB]. Both lead to the same results. 

For the liquid deuterium we fitted the total cross section vs. k using: 

A , , 

Otot = ^0 + (4-3) 

where = 27r/ A is the neutron wave vector in the material and (Jq can be interpreted 
as an energy independent incoherent elastic cross section. The total cross section 
Gtot is the sum of energy independent incoherent elastic cross section (Jq (4.1 b) and 
the capture plus inelastic cross section (up-scattering) which both are inversely pro- 
portional to the neutron wave vector. In Fig. 14. 6t the experimental data is plotted 
together with the fitted function. One notices that there is no visible deviation from 
the 1/f law. The parameter A which is a measure of inelastic cross section ys, A = 
(0.742 ± 0.009) b A" . The incoherent cross section value obtained from the fit is 
(To = (4.19 ± 0.11) b and agrees with the theoretical value to within one standard 
deviation. 
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Figure 4.5: The total cross section per deuterium molecule for liquid deuterium and 
solid deuterium plotted versus the neutron wave vector k. 
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Figure 4.6: The total cross section per deuterium molecule for liquid deuterium at 
19 K plotted versus the neutron wave vector k. The solid curve shows the result of a 
fit of the form atot = ctq + A/k. 
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The fitted function used for liquid deuterium is not suitable for the transmission 
data taken with the solid deuterium. Different dynamics of neutron scattering from 
the solid phase requires more advanced models, e.g. a model that includes the elastic 
coherent scattering due to inhomogeneities in the crystal. The fitted function which 
was applied to the data from the measurements of solid deuterium at 18 K follows the 
idea of Engelmann and Steyerl. In their paper they show that very cold neutron 
transmission can be used to investigate inhomogeneities in the matter f section [2.1.11) . 
The developed technique is alternative to the more conventional small angle scatter- 
ing. In this model the elastic scattering is described by the Eq. 12.201 of chapter 2. 
The total cross section may be written as: 



where g{kR6i) is given by Eq. 12.211 6'i = 20 ° is the detector aperture and R refers 
to the radius of submicroscopic structures. 



o-toi = do + — + 5 



g{kROi) 



(4.4) 
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Figure 4.7: The total cross section per deuterium molecule for solid deuterium at 
18 K plotted versus the neutron wave vector. The curve through the experimental 
points is a fit according to Eq. 14.41 
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Figure mZI shows the experimental data and the obtained fit. The incoherent elastic 
contribution do was taken from the measurement with liquid deuterium and A, B and 
R were varied. 

The parameter A = (0.266 ± 0.004) b A obtained from the fit is 2.8 times 
lower than that of liquid deuterium at a temperature of 19 K. This indicates strongly 
suppressed up-scattering for the solid deuterium as compared with liquid deuterium 
and agrees with a calculation presented by Keinert [32]: for the neutron energy of 
200 neV Keinert calculated the up-scattering rate for solid ortho-deuterium at tem- 
perature of 17 K to be 3.1 times lower than that of liquid deuterium at a temperature 
of 19 K. 

The interpretation of parameter B is the following 

B = M%''R^6p' (4.5) 

where h is the bound coherent scattering length, 6 p = p — po is the density difference 
to the surrounding medium. 

Parameter R refers to the radius of spherical inhomogeneities and the value ob- 
tained from the fit is = (9.6 ± 1.9) A. 

The same fitting procedure was applied to the data from the transmission mea- 
surements of the crystal at temperatures of 16 K, 14 K and 12 K. The experimental 
data and obtained fits are shown in Figs. 14. 8( 14.91 and I4.10( respectively. 
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Figure 4.8: The total cross section per deuterium molecule for SD2 at 16 K. 
curve through the experimental points is a fit according to Eq. 14.41 
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Fit parameters: 
chi'/ndf = 19.8/22 
A = (0.1686 ± 0.0040) b A"' 
B = (0.0014 ± 0.0003) b A"' 
R = (5.6 ± 1.1) A 
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Figure 4.9: The total cross section per deuterium molecule for SD2 at 14 K. The 
curve through the experimental points is a fit according to Eq. 14.41 
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B = (0.0005 ± 0.0003) b A"' 
R = (5.2± 1.3)A 
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Figure 4.10: The total cross section per deuterium molecule for SD2 at 12 K. The 
curve through the experimental points is a fit according to Eq. 14.41 
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The interpretation of the total cross section measurement for the deuterium crystal 
at temperatures below 12 K is also based on the model which includes both incoherent 
inelastic scattering (proportional to l/fc) and coherent elastic scattering (proportional 
to l/k"^)- The function which fits the best the data taken with the crystal at the 
temperatures below 12 K has a form: 

A B , 
a = ao + - + ^ (4.6) 

where parameter A is a measure of inelastic scattering and B gives an information 
about the strength of coherent elastic scattering. Figure 1 II shows that the exper- 
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Figure 4.11: The total cross section per deuterium molecule for solid deuterium at 
5 K vs. the neutron wave vector k. The curve through the experimental points is a 
fit according to Eq. 14.61 



imental point are well represented by the curve of the form 14.61 

The value for A = (0.056 ± 0.011) b A" is now about 3 times smaller than that 

obtained for crystal measurements at 18 K. 

The same analysis was applied to the data from the transmission measurements for 

crystal at temperatures of 10 K and 6 K. Figures [4 . 1 21 and l4 . 1 31 show the experimental 

data and the obtained fits for 10 K and 6 K, respectively. 
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Fit parameters: 
Chi' / ndf = 30.8 / 22 
A = (0.0635 ± 0.0044) b A"' 
B = (0.0032 ± 0.0003) b A"' 
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Figure 4.12: The total cross section per deuterium molecule for solid deuterium at 
10 K plotted versus the neutron wave vector. The curve through the experimental 
points is a fit according to Eq. 14.61 



We compared the experimental values of parameter A with the one phonon up- 
scattering cross sections (Eq. I2.38|l calculated for the initial neutron wavelength of 
60 A (or wave vector k = 0.1 A ). In Fig. 14.141 the experimental and theoretical 
values, given in b m s~^, are plotted versus the temperature of solid deuterium. It is 
seen that the experimental results reveal a different trend then theory. They agree 
well only for the crystal temperature of 18 K. 

To check the reproducibility of the measurement and sample preparation proce- 
dures we repeated the transmission measurements with another deuterium crystal 
with high ortho-deuterium concentration prepared in the same way. In Figs. 14.151 
and 14.161 the curves are plotted for both crystal samples at temperatures 18 K and 
5 K, respectively. We conclude that both measurements are consistent. 
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Figure 4.13: The total cross section per deuterium molecule for solid deuterium at 
6 K plotted versus the neutron wave vector. The curve through the experimental 
points is a fit according to Eq. 14.61 
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Figure 4.14: The up-scattering cross section vs. temperature for solid D2 as measured 
(red circles) and as follows from the one phonon up-scattering calculations (blue line). 
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Figure 4.15: The total cross section per deuterium molecule for solid deuterium at 
18 K for two different crystals prepared in the same way. 
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Figure 4.16: The total cross section per deuterium molecule for solid deuterium at 
5 K as obtained for two crystals prepared in the same way. 
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4.2.2 UCN data 

Figure 14.171 shows plots of the neutron transmission loss cross section versus neutron 
wave vector k for the liquid ortho-deuterium (cq = (98.2 ± 0.9)%) and solid ortho- 
deuterium at temperatures of 18 K and 5 K. The up-scattering cross section for the 
neutron wavelength of 550 A (or an energy of 270 neV) is 2.7 times smaller for solid 
deuterium at 18 K as compared with that for liquid deuterium and agrees with the 
value obtained from the VCN analysis. In Fig. 14.171 one also sees that, over the full 
energy range covered, at 5 K the inelastic cross section is a factor 2.2 smaller than 
at 18 K. This is 2 times less than the result obtained from the VCN analysis. 

Figure 14.181 shows the cross section as a function of the deuterium crystal tem- 
perature. One can notice the tendency of suppression of the total cross section with 
decreasing temperature of the crystal. However, the difference in the cross section 
values between 5 K and 10 K is not visible. It is noteworthy that the measurement of 
the UCN yield from the model of solid-deuterium UCN source performed by Serebrov 
et al. 1^ does not show any difference between the crystal at 10 K and 5 K either. 
Figure 14.191 shows the total cross sections for solid deuterium at temperatures of 
10 K and 5 K from both VCN and UCN measurements. It is seen that at the shorter 
wavelengths (30 A - 100 A) the cross section values for crystal at 10 K are higher 
as comapred with 5 K data, which may indicate that the coherent scattering plays a 
crucial role at the larger wavelengths. This effect is especially important at temper- 
atures below 10 K where the incoherent up-scattering cross section is relatively low. 
The cross section curves for solid deuterium at 18 K and 5 K are plotted in Fig. 14.201 
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Figure 4.17: The total cross section per deuterium molecule for liquid deuterium at 
a temperature of 19 K and solid deuterium at temperatures 18 K and 5 K. 
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Figure 4.18: The total cross section per deuterium molecule for solid deuterium at 
temperatures 18 K, 16 K, 14 K, 10 K and 5 K plotted vs. wave vector. 
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Figure 4.19: The total cross section per deuterium molecule for solid deuterium at 
temperatures 10 K and 5 Kas obtained from the measurements with VCN and UCN. 
In the wave vector range between 0.03 and 0.04 A the statistics is low thus the 
systematics effects dominate. This is the reason of increase of the cross section in 
this region. 
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Figure 4.20: The total cross section per deuterium molecule for solid deuterium at 
temperatures 18 K and 5 Kas obtained from the measurements with VCN and UCN. 



Chapter 5 
Conclusions 



For the construction of a high density superthermal UCN source several parameters 
need to be optimized. The production of UCN in the converter medium and the losses 
of UCN in this medium are defining the overall efficiency at which such a source can 
operate. Also the extraction efficiency from the source into an experiment itself is of 
great importance. 

To investigate the losses of UCN before they can exit the solid deuterium con- 
verter, neutron transmission measurements through solid deuterium samples have 
been performed with VCN and UCN beams. The storage lifetime related to these 
UCN losses in solid deuterium is dominated by four contributions PTJ: 

• the nuclear absorption of neutrons on deuterons, 

• temperature dependent phonon - UCN interaction in the deuterium crystal, 

• UCN up-scattering processes by contaminants such as H2, 

• the presence of para-D2 which behaves as a contaminant. 

The work presented in this thesis was devoted to the understanding of the tem- 
perature dependent losses namely losses due to up-scattering of UCN by picking up 
energy from phonons present in the deuterium. 

We have shown that the inelastic cross section for neutron wavelength range be- 
tween 40 A and 500 A in deuterium decreases by a factor 2.8 at the solidification of 
liquid deuterium f section 14.2.11) . This is in agreement with the calculations done by 
Keinert [HHl for a neutron wavelength of 640 A as well as with the results presented 
by Serebrov 37 . The experimental value of incoherent cross section obtained from 
the liquid deuterium data is (4.19 ± 0.11) b and agrees with the theoretical value to 
within one standard deviation. 

For solid deuterium the one phonon up-scattering is believed to be the only mech- 
anism which is responsible for temperature dependent losses of UCN (inelastic cross 
section). The interpretation method of the solid deuterium data that we used is 
based on a model developed by Steyerl and collaborators. It includes thermal inelas- 
tic scattering and the coherent elastic scattering on spherical inhomogeneities which 
are present in the material. Our measurements carried out with very cold neutrons 
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(in the wavelength range of 40 A - 190 A) give results that are consistent with this 
model. However, our data as obtained with ultra-cold neutrons (wavelength range 
210 A - 500 A) cannot be properly described by this model. The most likely expla- 
nation for this difference lies in the wavelength dependence of the neutron scattering 
on the material structure: the sensitivity to the scattering on inhomogeneities in the 
material is increased significantly at the larger wavelengths jHE]- 

From the fits according to the model mentioned above we were able to distinguish 
the inelastic part (parameter A) of the total cross section (section 14.2. 1|) for the tem- 
peratures of the crystal between 18 and 5 K. The comparison between up-scattering 
cross sections at different crystal temperatures shows that for wavelength range of 
40 A - 190 A: 

• At the temperature of 10 K the up-scattering cross section decreases 12 times as 
compared with liquid deuterium and 4.2 times as compared with solid deuterium 
at 18 K. 

• Further cooling of the crystal to the temperature of 5 K suppress the up- 
scattering only by 10%. 

At the wavelengths of 210 A - 500 A we observed that: 

• The up-scattering cross section for a crystal temperature of 10 K decreases 6 
times as compared with that for liquid deuterium and 2.2 times as compared 
with solid deuterium at 18 K. This is about 2 times less than the result obtained 
from the VCN analysis. 

• There is no visible reduction of inelastic cross section between measurements 
with the crystal at 10 K and 5 K. 

A possible explanation of those effects could be connected with a strongly enhanced 
coherent scattering of UCN. For a full understanding of this phenomenon more ex- 
tensive models are required. 

The results from neutron transmission studies are in accordance with optical in- 
vestigation of deuterium crystal and with Raman spectroscopy. From the pictures 
that were taken with the digital camera we can see that the transparency of crystal is 
changing with the temperature and the number of cracks is increasing with decreas- 
ing temperature. The analysis of Raman spectra shows that the orientation of the 
crystallites is changing during cooling which may influence the coherent scattering 
of neutrons. Clearly, measurements of this type will prove to be very important in 
further optimising the way SD2 is prepared and kept in the UCN source. 



Appendix A 



Time of flight spectrum and 
calibration 

A.l VCN experiment 

A. 1.1 Dead time correction of time of flight spectra 

The time-of-flight (TOF) spectra were measured for VCN packages that passed the 
neutron chopper. The raw spectra have to be corrected for the effect of dead time 
- specific for detectors and subsequent electronics. When the detecting system pro- 
cesses an event no additional events can be registered during a time interval of ap- 
proximately 20 (J.S. The influence of the detector dead time on the count rate can be 
calculated by analysing the distribution of time interval between subsequent events. 

The expected exponential distribution for an ideal (zero dead time) system follows 
from the Poisson statistics. The TOF spectra are divided into 2 ms bins and for each 
bin we build a distribution of time difference between subsequent events. Figure \Ka\ 
shows an example of this distribution for TOF between 20 and 22 ms taken in one 
run. 

The analysis consist of two steps. In the first one we determine the average dead 
time correction for every 2 ms bin: 

• An exponential function is fitted to an exponential function to the data between 
0.05 and 0.35 ms (Fig. Ell). 

• The integral of the fitted exponential distribution function extrapolated to zero 
time difference is calculated. 

• The ratio of the integral of the fitted function to the observed number of events 
is determined and it is this number that gives the dead time correction. 

In the second step the dead time loss factor as a function of TOF is interpolated 
using cubic splines (see Fig. IA.2|) . The interpolation is restricted to the values for 
which 0.6 < /Ndf < 1.7. l^df values refer to the fitted exponential functions 
from the first step. Figure IA.3I shows the effect of the dead time correction. The 
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Figure A.l: Time difference distribution of subsequent events for VCN's with time-of- 
flight between 20 and 22 ms with an exponential distribution curve y = Pq- exp (Pix) . 



dead time correction changes the shape of the spectrum and can reach 25% in the 
maximum of the TOF distribution. 



A.l. VCN EXPERIMENT 



57 



0.30 



0.24 



o 

-I—' 

^ 0.18 

05 
U3 

_o 

.i 0.12 
u 

CO 

o 
Q 

0.06 



0.00 



T ' 1 ' 1 ' I 



□ experimental data 
spline interpolation 




J 1 I 1 I 1 L 



10 15 20 25 30 35 40 45 
Time of fligth [ms] 



Figure A. 2: Dead time loss factor interpolated with cubic spline. 
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Figure A. 3: VCN TOF spectra without (blue solid line) and with (red dashed line) 
dead time correction for the empty cell. The bin width is 80 [is. 
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A. 1.2 Deconvolution of the time-of- flight spectra 

The raw TOF spectra are convolved with the chopper opening function as described in 
section IH. 1.61 To obtain the true TOF spectra a deconvolution procedure is necessary. 
We tried two methods (i) direct, based on FFT technique and (ii) iterative. The 
stable result in the direct deconvolution is strongly dependent on the input data 
as FFT techniques are sensitive to noise. Since the measured TOF spectra carry a 
non-negligible fraction of noise, the iteration method seems more appropriate. The 
iteration method assumes an analytical form of the unfolded function (i.e. the true 
spectrum) with free parameter and performs a convolution with the known chopper 
opening function. The resulting spectrum should be a good approximation to the raw 
spectrum as recorded by DAQ system. To optimize the fit, the following expression 
for chi-square is evaluated: 

X' = if{^)-yi^)f*wi^) (A.l) 

where f{i) is the folded theoretical function, y{i) is the value experimentally measured 
and w{i) is the weight of the data point (inverse square of the statistical error). The 
deconvolution factor is given by 

where g{i), f{i) are the values before and after the convolution respectively. The 
deconvolved spectra (with their natural statistical fluctuations) are then obtained by 
multiplying the raw data by the deconvolution factor. 

The analytic function used to parameterise the VCN TOF spectra is a sum of two 
modified Maxwellian time distributions and a constant background. 

f{x) =pi+ , ,„ exp ( -- — ^^^1+7 — Ml exp ( -- — ^^^-^ ) (A.3) 



{x-p2y V {'X-P2YJ {x-p2y^ V {'X-P2) 

Figure I A. 41 shows the experimental data and the fitted function /. The effect of 
deconvolution is presented in Fig. lA.51 The deconvolution factor d is shown in Fig. lA.6l 
and the deconvolved experimental spectrum is presented in Fig. I A. 71 
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Figure A. 4: An experimental time of flight spectrum fitted with the analytical func- 
tion /(red line). The green circles are the fit residuals. 
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Figure A. 5: Function / before deconvolution (red line) and after deconvolution g 
(green line). 
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Figure A. 6: A plot of the deconvolution factor g/f. 
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Figure A. 7: Time of flight spectra before and after deconvolution. 
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A. 1.3 Zero time calibration method 

The TOF spectrum needs to be corrected for an effect that arises because of a differ- 
ence in the time at which the chopper start signal is generated and the actual time 
at which the neutrons start passing through the chopper. This effect will appear 
as a shift in the TOF spectrum. This section describes the calibration method to 
determine the zero point on the time scale. 

For calibration, data was collected with the detector placed at two positions: 
"front" (F) and "extended" (E), with the same coUimation (section 13.1.11) . The 
difference in the flight path between these two configurations is sef = 511 mm. 
Figure IA.8I shows the time-of-fiight spectra recorded at these two positions. 




Time of flight [ ms ] 



Figure A. 8: Time of flight spectra for "front" and "extended" detector positions. 
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Figure A. 9: A beam profile as measured for tlie "front" detector position. 




Figure A. 10: A beam profile as measured for the "extended" beam position. 
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Assuming that the mean arrival time does not change with spreading of the TOF 
spectrum, the mean velocity can be calculated by dividing the distance sef by the 
difference in the mean arrival time at both positions: 



ts — tp 



(A.4) 



A true mean arrival time at the "front" position is then determined by dividing the 
length of the flight path [sp = 1436 mm ) by the mean velocity: 



V 



(A.5) 



The correction to the arrival time was than calculated by subtracting the measured 
time from the true time. 

There are two contributions to the systematic error connected with this calibration 
method: 

• gravitational and beam divergence effects, 

• determination of the distance with the assumption that neutrons react with the 
^He in the middle of the detector. 



Collimator 



I 



Sample Detector 

Figure A. 11: A schema showing the dimensions needed for calculations of the beam 
spot. 



To minimise this uncertainty, the TOF spectra have been built for the events falling 
into the collimator projected spot area defined by the spot size parameter b: 

b = d*(^^ + ?j (A.6) 

where d = 20 mm is the diameter of the collimator, / = 1260 mm is the collimation 
length and x is the distance from the collimator to the detector (see Fig. lA.llI) . Thus 
the beam spot size is 26 mm for the "front" detector localization (x = 174 mm) and 
42 mm for the "extended" detector position {x = 685 mm). Figures fA.9l and lA.lOl 
show the beam profiles for the "front" and "extended" detector position. Divergence 
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of the beam and gravitation cause missing of events at the "extended" detector po- 
sition, and consequently, changing of the mean arrival time. The contribution from 
this effect to the calibration error was estimated to be of 0.04 ms. The uncertainty 
in determining the exact flight path length gives a systematic error of 0.28 ms. 

The obtained correction to the arrival time (shift in the TOF spectrum) is (0.80 
± 0.28) ms. The effect of calibration is shown in Fig. IA.121 one curve shows a TOF 
spectrum before and the other one after calibration. 

In order to check the reliability of this method we calculate the time-of-flight distri- 




Time of flight [ ms ] 

Figure A. 12: Time of flight spectra before and after calibration. 

bution at the close position from known distribution at the "extended" position. The 
comparison is shown in Fig. IA.13I 

A. 1.4 Background measurements 

The background rate for the VCN experiment was measured to be about 2 events 
per second. The measurement was performed with the main shutter closed when the 
reactor was running. The main contribution comes directly from the turbine. 
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Figure A. 13: Comparison between the measured spectrum at the "front" position 
(sohd hne) and transformed to that position from a measurement at the "extended" 
position (dashed hne). 
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A. 2 UCN experiment 

A. 2.1 Dead time consideration 

The neutron flux at the UCN test position is considerably lower than that at the 
VCN beam: the UCN flux of the order of 20 cm~^ s~^ for v = 13 m/s and about 
4000 cm~^ s~^ at the velocity of 70 m/s. Consequently, for the experiments carried 
out with the UCN beam, no dead time corrections are required. 



A. 2. 2 Calibration of the timing in the TOF spectra 

As is the case of the VCN experimental setup, the neutron chopper operating in the 
UCN setup generates a start signal that does not exactly coincide with the opening 
of the neutron window itself f section l3.1.6p . This causes shift in the TOF spectra and 
needs to be corrected by a calibration procedure. Two independent methods have 
been developed to calibrate UCN time-of-flight spectra. 

The first one is based on a direct experimental determination of the number of 
steps the stepping motor moves between the moment the TTL start signal appears 
and the moment the neutron start passing through the chopper. The chopper is driven 
slowly, step by step, by a control program running in the Lab View environment. 
This is why the exact angular position of the stepping motor at which the TTL 
signal appears can be determined. Likewise by driving the chopper further, step 
by step, the exact angular position of the stepping motor at which the UCN pass 
through the chopper and are detected by the detector is determined, too. One thus 
knows the exact number of steps Nsw between the TTL start signal and the chopper 
opening. Knowing also the frequency u at which the stepping motor is driven one 
can then calculate the time offset to to be applied to raw TOF spectra. In a similar 
procedure, one can also experimentally determine the number of steps Nqc of the 
motor between opening and closing the chopper window as well as the number of 
steps of a full chopper cycle. This method was applied for the measurements with 
the beam collimator and without it. The obtained results are collected in Tab. lA.ll 



Setup 


Nsw 


Noc 


Full chopper turn 


With collimator 


591 


60 


864 


Without collimator 


595 


52 


864 



Table A.l: Results from the calibration measurements 



In the calculations we have assumed that t = in the centre of the chopper 
opening. For the measurements with the chopper speed of 2500 steps/s (3 Hz) the 
calculated offset is: 

to = = 248 ms (A.7) 

The precision of this results, limited by the uncertainty in counting the number of 
steps performed, is 1 ms. 
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The other method of cahbrating the TOF spectra is based on measuring the arrival 
time of the fastest (slowest) neutrons at various chopper speeds and fitting a linear 
function to the data. Examples for the measurements without collimator are plotted 
in Figs IA.14I and IA.15I The time offset between TTL signal and chopper window 
opening is then given by 



where Pi is the fit parameter and u is the chopper frequency. The calculated offset 
(average from both measurements) is to = (250 ± 2) ms. 

Both these calibration methods rely on a precise knowledge of the chopper speed. 
The chopper has been tested and has shown to be operating stable over long periods. 
The chopper speed scales linearly with the true disks speed. Both calibration methods 
gave consistent results. Fi gures [A . 1 61 and I A . 1 71 show a time-of- flight spectrum with no 
calibration applied and a plot of the same spectrum with proper timing calibration. 
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Figure A. 14: The calibration of the UCN TOF spectrum based on the arrival time 
of the slowest neutron. 
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Figure A. 15: The calibration of the UCN TOF spectrum based on the arrival time 
of the fastest neutron. 



A. 2. 3 Background measurements 

The background for the UCN experiment appears to be of the order of 0.1 counts per 
second. It was measured with the main shutter closed and the reactor running. 
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Figure A. 16: A time-of- flight spectrum before calibration. The spectrum is "split" 
into two parts by the start signal which is generated during the chopper window is 
open. The bin width is 4 ms. 
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APPENDIX A. TIME OF FLIGHT SPECTRUM AND CALIBRATION 
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Figure A. 17: A time-of-flight spectrum with time calibration. The bin width is 4 ms. 
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